Phase growth in bistable systems with impurities 
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A system of coupled chaotic bistable maps on a lattice with randomly distributed impurities is 
investigated as a model for studying the phenomenon of phase growth in nonuniform media. The 
statistical properties of the system are characterized by means of the average size of spatial domains 
of equivalent spin variables that define the phases. It is found that the rate at which phase domains 
grow becomes smaller when impurities are present and that the average size of the resulting domains 
in the inhomogeneous state of the system decreases when the density of impurities is increased. The 
phase diagram showing regions where homogeneous, heterogeneous, and chessboard patterns occur 
on the space of parameters of the system is obtained. A critical boundary that separates the regime 
of slow growth of domains from the regime of fast growth in the heterogeneous region of the phase 
diagram is calculated. The transition between these two growth regimes is explained in terms of 
the stability properties of the local phase configurations. Our results show that the inclusion of 
spatial inhomogeneities can be used as a control mechanism for the size and growth velocity of 
phase domains forming in spatiotemporal systems. 

PACS numbers: 05.45.-a, 89.75.Kd 



I. INTRODUCTION 

There has been much recent interest in the study of 
spatiotemporal dynamical processes on nonuniform or 
complex networks. In this context, coupled map lattices 
[l| have provided fruitful and computationally efficient 
models for the investigation of a variety of dynamical pro- 
cesses in spatially distributed systems. In particular, the 
discrete-space character of coupled map systems makes 
them specially appropriate for the investigation of spa- 
tiotemporal dynamics on nonuniform networks that can 
represent models of heterogeneous media. The nonuni- 
formity may be due to the intrinsic heterogeneous nature 
of the substratum, typical of pattern formation in biolog- 
ical contexts, or it may arise from random imperfections 
or fluctuations in the medium at some length scales. Such 
heterogeneities can have significant effects on the forma- 
tion of spatial patterns, for example, they can induce 
reberverators in excitable media and defects can serve as 
nucleation sites for domain growth processes. 

Recently, the study of the phase-ordering properties of 
systems of coupled chaotic maps and their relationship 
with Ising models in statistical physics has been a focus 
of attention [3, H, HIS S 0, Hi • These works have mainly 
assumed the phase competition dynamics taking place on 
an ordered spatial support. This article investigates the 
phenomenon of phase growth in coupled chaotic maps on 
lattices with randomly distributed impurities as a model 
for studying this process on nonuniform media. In partic- 
ular, this model of coupled maps on nonuniform networks 
yields a scenario to explore the role that the local config- 
urations of the underlying lattice play on the statistical 
properties of phase ordering processes on spatiotemporal 
systems. In addition, this model can be seen as a sim- 



pler but computationally more efficient alternative to the 
study of phase separation phenomena on media with im- 
purities than conventional computational fluid dynamics 
techniques. 

In Sec. II, we present a procedure for the construc- 
tion of a lattice with impurities and define the coupled 
map model on this network. The phase growth dynamics 
in the presence of impurities is studied in Sec. III. The 
phase diagram of the system on the space of its parame- 
ters is obtained in this Section. Section IV contains the 
conclusions of this work. 



II. COUPLED MAP LATTICE WITH 
IMPURITIES 



To generate a lattice with impurities, we start from a 
two-dimensional array of cells of size L x L with periodic 
boundary conditions and remove a given fraction p of 
cells at random. The removed sites can be considered 
as impurities, defects, or random imperfections on the 
spatial support. We define an impurity as a non- active 
site, i.e., a site that possesses no dynamics. The fraction 
of impurities in the lattice can be characterized in terms 
of the minimum Euclidean distance d between impurities, 
as shown in Figure 1. In this way, the density p scales 
with the minimum distance between impurities as p ^ 
0.625(i~^. We shall use values d > 2 in order to avoid 
contiguous impurities on the lattice. 

The N = L^{1 — p) active sites can be enumerated 
by an index i — I, . . . , N. The equations describing the 
dynamics of the diffusively coupled map system defined 
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FIG. 1: (Left) Spatial support showing active sites (o) and 
impurities (•). Impurities are randomly placed keeping a 
minimum distance d between them. (Right) The density of 
impurities scales as p = 0.625d~^. 

on such nonuniform lattice are 

:rJ+i = (l-6)/(xj) + -^^/(a;^), (1) 

where gives the state of an active site i at time t; 
v"^ is the set of the nearest active neighbors (von Neu- 
mann neighborhood) of site z, and A/"' G {1, 2, 3, 4} is the 
cardinality of v"^] the parameter e G [0, 1] measures the 
coupling strength, and f{x) is a map that expresses the 
local dynamics. Impurities are not coupled to any other 
cell in the lattice. 

The above coupled map equations on a lattice with 
impurities can be generalized to include other cou- 
pling schemes and neighborhoods, higher dimensions or 
continuous-time local dynamics. Different spatiotcmpo- 
ral phenomena can be also be studied on such nonuniform 
structures by providing appropriate local dynamics and 
couplings. 

In order to describe a bistable dynamics, we assume a 
piecewise local map Q 

( -2^/3 -^ix, if a; e [-1,-1/3] 
fix) = \ ^ix, if .TG [-1/3,1/3] (2) 
[ 2^/3 -Aia;, if x e [1/3,1] . 

When the parameter fi £ (1,2) the map has two symmet- 
ric chaotic band attractors, one with values xl > and 
the other with xl < 0, separated by a finite gap about the 
origen. Then the local states have two well defined sym- 
metric phases that can be characterized by spin variables 
defined as the sign of the state at time t, al — s\g\\{x\). 

III. PHASE GROWTH IN THE PRESENCE OF 
IMPURITIES 

We fix the local parameter at /i = 1.9 and set the ini- 
tial conditions as follows: one half of the active sites are 
randomly chosen and assigned random values uniformly 
distributed on the positive attractor while the other half 
are similarly assigned values on the negative attractor. If 
the number of active sites is odd, then the state of the 
remaining site is assigned at random on either attractor. 



In regular lattices (p = 0) phase growth occurs for 
values e > e,,, where eo — 0.67 |2]. In contrast, in the 
medium with impurities there exist a minimum value of 
p for which the domains formed by the two phases reach 
a frozen configuration for all values of the coupling e. 
Figure 2 shows stationary patterns emerging in the sys- 
tem for different values of parameters. The top panels 
show that the average size of domains decreases when 
the density of impurities in the system is increased. The 
bottom panels reveal the presence of complex domains 
where both phases coexists in a chessboard (also called 
antiferromagnetic) configuration for large values of the 
coupling strength. 




FIG. 2; Stationary patterns on a lattice of size 100 x 100 
with impurities, for different values of parameters . Top left: 
p = 0.0097 (d = 8), e = 0.73. Top right: p = 0.039 (d = 4), 
e = 0.73. Bottom left: p = 0.0097 (d = 8), e = 0.95. Bottom 
right: p = 0.0097 (d = 8), e = 0.98. 

To characterize the phase ordering properties of the 
system Eq. ([T|) we use the normalized size of one phase 
domain as a function of time as an order parameter, de- 
fined as 

L/2 

r— 1 i^j 

where is the Euclidean distance between sites i and j. 
Figure 3 shows the average asymptotic value {Roo) as a 
function of e in the absence of impurities. The quantity 
{Rcx)) exhibits a continuous phase transition from a het- 
erogeneous state, characterized by a small value of (i?oo), 
to a homogeneous state, for which (Roo) = 1, at the crit- 
ical value Eo of the coupling parameter reported in [2]. 

Before reaching its stationary value, the average do- 
main size increases in time as (i?t) t", where the expo- 
nent a characterizes the rate of phase growth. Figure 4 
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FIG. 3: Asymptotic (R) as a function of e for d = cxd (p = 0), 
averaged over 40 realizations of initial conditions. Error bars 
correspond to the standard deviations. 



shows the evolution of (Rt) with fixed value of e for dif- 
ferent values of the distance d between impurities. We 
observe that the phase growth process follows two well 
differentiated regimes during its time evolution. At early 
times t ^ 100, we find {Rt) ^ with an exponent 
ai w 0.071 whose value is not appreciably affected by 
the presence of impurities. For later times, the scaling 
behavior changes to (Rt) t"^, with a2 > ai indicating 
that domains grow faster in this regime, and where a2 de- 
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FIG. 4: Log-log plot of the average domain size {Rt} vs. t for 
d = CO (p = 0) (empty circles), and d = 6 (p = 0.017) (solid 
circles), with fixed e = 0.73, averaged over 40 realizations. 
The values of the scaling exponents ai and a2, as well as a 
typical error bar are indicated on each curve. 



pends on the density of impurities. Impurities inhibit the 
growth of the phase domains and prevent the system from 
reaching a homogeneous state for values of e > Eq. Fig- 
ure 4 shows that the average size of the resulting domains 
in the inhomogeneous state are smaller in the presence 
of impurities, in agreement with the behavior observed 
in other bistable systems on media with impurities Tol. 
Theoretical models and experiments with binary fluids 
have also shown that phase growth is restricted when 
impurities are present [ll|, [I4I 

The domain growth regimes characterized by the ex- 
ponents ai and 012 in Fig. 4 can be understood in terms 
of the stability of local configurations of the two phases. 
With this aim, we define the fraction of sites in a given 
phase that have k neighbors in that same phase at time 
t, given by 



1 ^ 



with fc = 0, 1,2,3,4; where 



0, if i:,eu^5^,..,^k 



(4) 



(5) 



Note that X]fe^t(^) — 1- Figure 5 shows the local 
spatial configurations corresponding io k = 0,1,2,3,4. 
The patterns associated to k — 0,1,2,3 represent config- 
urations where the two phases are in contact, while the 
homogeneous state corresponds io k ~ A. 
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FIG. 5: Local configurations taken into account by the frac- 
tions G{k), for fc = 0, 1, 2, 3, and 4. 

Figure 6 shows the time evolution of the fractions 
Gt{k) for two different values of d. We observe that in 
both cases Gt{0) and G't(l) vanish at t ~ 100. Thus, the 
decay of in the number of local configurations associated 
to A: = and k = 1 determines the regime of slow do- 
main growth characterized by the exponent ai in Fig. 4. 
Phase growth in this earlier regime occurs by the addi- 
tion of single cells to domains whose size is of the order 
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FIG. 7: Top: Mechanism of fast phase growth when two do- 
mains of the same phase join together to form a larger domain. 
The arrow indicates the direction of time. Bottom: The local 
configuration changes occurring in the boxes marked in the 
top panel are illustrated. 
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FIG. 6: Fractions Gt{k) vs. t for different values of d, with 
fixed e — 0.73. Gt{0) (continuous line); G(l) (long-dashed 
line); Gt(2) (short-dashed line); and Gt(3) (dotted line), (a) 
d — oo. (b) d = 6. 



of one cell. The spatial scale involved in this process is 
smaller than the minimum distance between impurities d 
and, therefore, the presence of impurities does not affect 
this mechanism of phase growth. As a consequence, the 
exponent ai does not depend on d. 

For t > 100, the fractions Gt{2) and Gt(3) keep de- 
caying until they reach their asymptotic values. Phase 
growth in this second regime occurs when domains in 
the same phase and whose average sizes are greater than 
one cell enter in contact. This process is illustrated in 
Fig. 7. At the local level, when two domains of equal 
phase get separated by a one-cell thick layer of cells in the 
opposite phase, there are cells on this layer that posses 
a local configuration described by A: = 2. The coupled 
map dynamics, Eqs. ([1]), acting on a cell in such local 
configuration will produce a change of its phase, yielding 



first another configuration of type k — 2, then transform- 
ing into a configuration of type k — S, and finally into 
the configuration associated to A; = 4. As a result, the 
domains join together forming a domain whose average 
size is much larger than the sizes of the initial domains. 
Through this growth mechanism, the average size of do- 
mains in the system increases faster than by the process 
the successive additions of single cells that takes place 
at earlier times. The difference in the growth velocity of 
phase domains in these two regimes is manifested by the 
fact that a2 > ai. 

When impurities are present, the average size of do- 
mains is comparable to the minimum distance d. Thus, 
there is a high probability that some impurities lie on 
the interface. These impurities constrain the transfor- 
mations of the local configurations k = 2 into k = 3, 
and into fc = 4 described above. Consequently, the phase 
growth process becomes slower, and the fractions Gt{2) 
and Gt(3) reach greater asymptotic values in the presence 
of impurities, as manifested in Fig. 6(b). Therefore, the 
asymptotic average size of domains for p is smaller 
than that for p = 0, as it was shown in Fig. 4. 

Figure 8 shows the growth exponents ai and a2 as 
functions of the coupling parameter e, for two different 
values of the density of impurities p. We observe that ai 
increases slowly with e and, as it was already manifested 
in Fig. 4, its behavior is unaffected by the presence of 
impurities. In contrast, the exponent a2 appears above 
some threshold value Cc > Eo and its behavior depends 
on both e and p. 

Figure 9 shows the growth exponent 012 calculated on 
the space of parameters (e, d). There is a critical bound- 
ary ec{d) that separates the regimes where fast growth 
occurs, characterized by a2 > 0, from the regime where 
only slow growth takes place, corresponding to a2 = 0, 
on the plane (e, d). Near the critical boundary edd), 
the exponent 012 can be described by the scaling relation 
a2 ~ (e — ec)'^, where 7 also depends on d or, equivalently, 
on the density of impurities p. 

Figure 10 shows the critical exponent 7 versus d. The 
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FIG. 8: Phase growth exponents ai (dashed Hne) and 0:2 (con- 
tinuous lines) vs. e for different values of d. The values shown 
are averages obtained over 40 realizations of initial conditions 
foe each value of e. Empty circles: ai for d — 00 {p = 0); 
empty squares: ai for d = 6 (p = 0.017); solid squares: 02 
for d = 00; solid circles: 02 for d = 6. Typical error bars are 
shown on each curve. 



exponent 7 can be well fitted by the scaling relation 
7(d) ~ (d - dcf, where 4 « 4.0 and 9 « 0.45. Be- 
low dc « 4.0 {p = 0.039) there is no regime of fast phase 
growth. 

Our results show that it is possible to control both the 
growth velocity and the size of phase domains in bistable 
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FIG. 10: The critical exponent 7 as a function of the distance 
between impurities d, obtained as averages over 40 realizar 
tions of initial conditions for each parameter value. Error 
bars are shown. 



media by the inclusion of impurities in the spatial sup- 
port of the system. When the density of impurities in- 
creases above some critical value for a given strength of 
the coupling, both the velocity at which domains grow 
and their stationary sizes reach smaller values than those 
corresponding to the absence of impurities. 

To characterize the emergence of chessboard domains 
for large values of the coupling e, as those shown in the 
bottom panels of Fig. 2, we calculate the asymptotic 
fractions Goo{k) as functions of e in Fig. 11. Note that 




FIG. 9: The fast phase growth exponent 02 as a function of 
e and d. The critical boundary ec(d) that separates the fast 
growth regime from slow growth regime is indicated by a thick 
continuous line on the plane (e, d). The values are obtained 
as in Fig. 8. 
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FIG. 11: Fractions Goo(fc) versus e for p = 0. Goo(O) (con- 
tinuous line) ; Goo(l) (long-dashed line); GoaCi) (short-dashed 
line); and Goo (3) (dotted line). Values shown correspond to 
averages over 40 realizations of initial conditions. 
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Goo(O), which vanished rapidly for small values of e, be- 
comes different from zero and increases for e > 0.9. The 
local configuration corresponding to fc = is effectively 
associated to a chessboard pattern. Thus the appear- 
ance of a finite value of the fraction Goo (0) at this large 
value of e signals the onset of chessboard domains in the 
system. Furthermore, we have found that the presence 
of impurities has little effect on the behavior of Goo(O), 
suggesting that the existence of chessboard domains is 
mainly associated to large values of the coupling param- 
eter and to the properties of the local dynamics. 

Figure 12 summarizes the collective behavior of the 
coupled map lattice with bistable dynamics in the pres- 
ence of impurities, Eq. H]). This figure shows the phase 
diagram of this system on the space of parameters (e, p). 
The regions where the homogeneous (one phase) and het- 
erogeneous (two coexisting phases), as well as the region 
where chessboard (antiferromagnetic) states occur are in- 
dicated. The critical boundary ec{p) for the onset of the 
regime for fast growth of domains (a2 > 0) on this plane 
is marked by a dashed line. 
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model to study phase ordering processes on media with 
impurities. The density of impurities p is related to 
the minimum distance between impurities d. We have 
shown that the rate at which phase domains grow be- 
comes smaller when impurities are present in the system. 
On the other hand, the average size of the resulting do- 
mains in the inhomogeneous state of the system decreases 
when the density of impurities is increased. 

We have calculated the critical boundary that sepa- 
rates the regime of slow growth of domains from the 
regime of fast growth on the space of parameters of the 
system (e,d). Along this critical boundary, the expo- 
nent for fast growth a.2 exhibits scaling properties along 
both parameter axes, manifested by the existence of the 
critical exponents 7 and B in Figs. 9 and 10. We have ex- 
plained the transition between these two growth regimes 
in terms of the stability properties of the local configu- 
ration measured by the fractions Gt(fc). 

For large values of the coupling strength e, the sys- 
tem displays a complex heterogeneous state consisting of 
domains of the two phases coexisting with spatial do- 
mains having a chessboard (antiferromagnetic) configu- 
ration. We have found that the appearance of chessboard 
patterns is not appreciably affected by the presence of im- 
purities. The phase diagram showing the regions where 
the different behaviors of the system occur on the space 
of parameters (e, p) was obtained. 

Our results indicate that the inclusion of impurities can 
be used as a control mechanism for the size and growth 
velocity of domains forming in bistable media. Recent 
studies of chaotic maps on complex networks [1, [l^ indi- 
cate that topology may play a decisive role in determining 
emerging collective behaviors. The present results sug- 
gest that spatial inhomogeneities may also be employed 
as a selection mechanism for patterns arising in general 
spatiotemporal systems. 



FIG. 12: Phase diagram of the system Eq. ([1} on the space 
of parameters (e, p). The regions where homogeneous, hetero- 
geneous, and chessboard phases exist are indicated by labels. 
The dashed line indicates the critical boundary for the on- 
set of the regime for fast growth of domains characterized by 
02 > 0. The critical value to for the occurrence of domain 
growth in the absence of impurities is marked on the vertical 
axis. 
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